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Abstract 

We propose a useful approach for investigating the statistical properties of foreign currency ex- 
change rates. Our approach is based on queueing theory, particularly, the so-called renewal-reward 
theorem. For the first passage processes of the Sony Bank US dollar / Japanese yen (USD/JPY) 
exchange rate, we evaluate the average waiting time which is defined as the average time that 
customers have to wait between any instant when they want to observe the rate (e.g. when they 
log in to their computer systems) and the next rate change. We find that the assumption of expo- 
nential distribution for the first-passage process should be rejected and that a Weibull distribution 
seems more suitable for explaining the stochastic process of the Sony Bank rate. Our approach also 
enables us to evaluate the expected reward for customers, i.e. one can predict how long customers 
must wait and how much reward they will obtain by the next price change after they log in to 
their computer systems. We check the validity of our prediction by comparing it with empirical 
data analysis. 
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I. INTRODUCTION 



Recently, internet trading has become very popular. Obviously, the rate (or price) change 
of the trading behaves according to some unknown stochastic processes, and numerous 
studies have been conducted to reveal the statistical properties of its nonlinear dynamics 
nj, y, yl. In fact, several authors have analysed tick-by-tick data of price changes including 

nnnpnnn 

the currency exchange rate in financial markets [4, l5|, la, 17|, l8|, l9|, Il0| . Some of these studies 
are restricted to the stochastic variables of price changes (returns) and most of them are 
specified by terms such as the fat or heavy tails of distributions l| . However, fluctuation 
in time intervals, namely, the duration in the point process ll|, ll2j might also contain 
important market information, and it is worthwhile to investigate these properties. 

Such fluctuations in the time intervals between events are not unique to price changes in 
financial markets but are also very common in the real world. In fact, it is wellknown that 
the spike train of a single neuron in the human brain is regarded as a time series, wherein 
the difference between two successive spikes is not constant but fluctuates. The stochastic 
process specified by the so-called inter-spike intervals (ISI) is one such example Q- 
The average ISI is of the order of a few milli-second and the distribution of the intervals is 
well- described by a Gamma distribution jl^ . 

On the other hand, in financial markets, for instance, the time interval between two 
consecutive transactions of Bund futures (Bund is the German word for bond) and BTP 
futures (BTPs are middle- and long-term Italian Government bonds with fixed interest 



rates) traded at London International Financial Futures and Options 



jxchange (LIFFE) is 

10 seconds and is well-fitted by the Mittag-Leffler function [6|, |7|, [ISj . The Mittag-Leffler 
function behaves as a stretched exponential distribution for short time-interval regimes, 
whereas for the long time-interval regimes, the function has a power-law tail. Thus, the 
behaviour of the distribution described by the Mittag-Leffler function changes from the 
stretched exponential to the power-law at some critical point [igI. However, it is nontrivial 
to determine whether the Mittag-Leffler function supports any other kind of market data, 
e.g. the market data filtered by a rate window. 

Similar to the stochastic processes of price change in financial markets, the US dol- 



lar/Japanese yen (USD/JPY) exchange rate of Sony Bank 



17|, which is an internet-based 



bank, reproduces its rate by using a rate window with a width of 0.1 yen for its individual 
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customers in Japan. That is, if the USD/JPY market rate changes by more than 0.1 yen, 
the USD/JPY Sony Bank rate is updated to the market rate. In this sense, it is possible for 



us to say that the procedure of determination o: 



essentially a first-passage process 
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;he Sony Bank USD/JPY exchange rate is 



22|. 



In this paper, we analyse the average time interval that a customer must wait until the 
next price (rate) change after they log in to their computer systems. Empirical data analysis 
has shown that the average time interval between rate changes is one of the most impor- 
tant statistics for understanding market behaviour. However, as internet trading becomes 
popular, customers would be more interested in the average waiting time — defined as the 
average time that customers have to wait between any instant and the next price change — 
when they want to observe the rate, e.g. when they log in to their computer systems rather 
than the average time interval between rate changes. To evaluate the average waiting time, 
we use the so-called renewal-reward theorem which is wellknown in the field of queueing 
theory [23|, . In addition, we provide a simple formula to evaluate the expected reward 
for customers. In particular, we investigate these important quantities for the Sony Bank 
USD/JPY exchange rate by analysing a simple probabilistic model and computer simula- 
tions that are based on the ARCH (autoregressive conditional heteroscedasticity) 25| and 



GARCH (generalised ARCH) Q, l26|, 27 1 stochastic models with the assistance of empirical 
data analysis of the Sony Bank rate X lo|. 

This paper is organised as follows. In the next section, we explain the method being 
used by the Sony Bank and introduce several studies concerning empirical data about the 
rate. In Sec. Ill, we introduce a general formula to evaluate average waiting time using the 
renewal-reward theorem and calculate it with regard to Sony Bank customers. Recently, one 
of the authors lo| provided evidence implying that the first-passage time (FPT) distribution 
of the Sony Bank rate obeys the Weihull distribution [28]. This conjecture is regarded as a 
counter part of studies that suggest that the FPT should follow an exponential distribution 
(see example in 291]). In the same section, we evaluate the average waiting time while 
assuming that FPT obeys an exponential distribution. Next, we com par e it with the result 
for the WeibuU distribution and perform empirical data analyses jo], [l^ . We find that the 
assumption of exponential distributions on the first-passage process should be rejected and 
that a WeibuU distribution seems more suitable for explaining the first-passage processes of 
the Sony Bank rate. Thus, we can predict how long customers wait and how many returns 
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they obtain until the next rate change after they log in to their computer systems. In Sec. 
IV, to investigate the effects of the rate window of the Sony Bank, we introduce the ARCH 
and GARCH models to reproduce the raw data before filtering it through the rate window. 
In Sec. V, we evaluate the reward that customers can expect to obtain after they log in to 
their computer systems. The last section provides a summary and discussion. 



II. DATA: THE SONY BANK USD/JPY EXCHANGE RATE 



The Sony Bank rate [l7j is the foreign exchange rate that the Sony Bank offers with 
reference to the market rate and not their customers' orders. In FIG. [H we show a typical 
update of the Sony Bank rate. If the USD/JPY market rate changes by greater than or 




FIG. 1: The Sony Bank USD/JPY exchange rate. The units of the horizontal and vertical axes are ticks 
and yen/dollar, respectively. 

equal to 0.1 yen, the Sony Bank USD/JPY rate is updated to the market rate. In FIG. 121 
we show the method of generating the Sony Bank rate from the market rate. In this sense, 
the Sony Bank rate can be regarded as a kind of first-passage processes. In TABLE [H we 
show data concerning the Sony Bank USD/JPY rate vs. tick-by-tick data for the USD/JPY 
rate from Bloomberg L.P. From these tables and figures, an important question might arise. 
Namely; how long should Sony Bank customers wait between observing the price and the 
next price change? This type of question never arises in the case of ISI or Bund futures, 
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FIG. 2: An illustration of generating the filtered rate (black circle) by the rate window with a constant 
width 2e in time (shaded area) from the market rate (solid line) . The units of the horizontal and the vertical 
axes are ticks and yen, respectively. 

TABLE I: The Sony Bank USD/JPY rate vs. tick-by-tick data for the USD/JPY rate 

Sony Bank rate tick-by-tick data 

Number of data per day ~ 70 ~ 10, 000 

The smallest price change 0.1 yen 0.01 yen 

Average interval between data ~ 20 minutes ~ 7 seconds 

because the average time intervals are too short to evaluate such measures. We would like 
to stress that in this paper, we do not discuss the market data underlying the Sony Bank 
rate; however, as we will see in the following sections, the main result obtained in this study, 
i.e. how long a trader must wait for the next price change of the Sony Bank rate from the 
time she or he logs on to the internet, is not affected by this lack of information. 

III. RENEWAL-REWARD THEOREM AND AVERAGE WAITING TIME 

From TABLE [H we find that the number of data per day is remarkably reduced from 
~ 10, 000 to ~ 70 because of an effect of the rate window of 0.1-yen width. As a result, the 
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average interval of exchange rate updates is extended to ~ 20 min. This quantity is one of 
the most important measures for the market, however, customers might seek information 
about the average waiting time, which is defined as the average time interval that customers 
have to wait until the next change of the Sony Bank USD/JPY rate after they log in to 
their computer systems. To evaluate the average waiting time, we use the renewal-reward 



theorem, which is wellknown in the field of queueing theory 23|, |2J] . We briefly explain the 
theorem below. 

Let us define N{t) as the number of rate changes within the interval (0, r] and suppose 
that a customer logs in to her/his computer system at time r {tn{t) < t < tn{t)+i)- Then, 
we defined the following quantity: 

W{r) = rjv(r)+i - T, (1) 

which denotes the waiting time for the customer until the next update of the Sony Bank 
USD/JPY rates when she or he logs in to a computer system at time r. Then, the renewal- 



reward theorem 



23 



24| implies that the average waiting time w can be written in terms of 
the first two moments, E{t) and E{t^), of the first-passage time distribution P{t) as 

where E{- ■ ■) denotes dt(- ■ ■)P{t). Therefore, if we obtain the explicit form of the FPT 
distribution P(t), we can evaluate the average waiting time w by using this theorem ([2]). 

The proof of the theorem is quite simple. Let us suppose that the points at which the 
Sony Bank USD/JPY rate changes are given by the time sequence ti,T2, - ■ ■ , tn{t)- For these 
data points, the first-passage time series is given by definition as to = ^1,^1 = T2 — ^1,^2 = 

— T2, - ' ' )^Af(r)-i = TAf(r) " tn(t)-i- Then, wc obscrvc that the time integral fQW{s)ds 
appearing in equation ([2]) is identical to ^f=Q\tiy /2, where (ti)^/2 corresponds to the area 
of a triangle with sides ti,ti and V2ti. As a result, we obtain 

wherein we used the fact that the expectation of the waiting time ti is given by 

1 1 r 



Therefore, 

1„„ ^ = ' (5) 

r-oo r E{t) ^ ' 

Thus, equation ([2]) holds true. 

Therefore, if we obtain the exphcit form of the first-passage time distribution the 
average waiting time w can be evaluated by theorem ([2]). However, for estimating the distri- 
bution we can only perform empirical data analysis of the Sony Bank rate. Apparently, 
the first-passage process of the Sony Bank rate depends on the stochastic processes of raw 
market data underlying the data that is available on the website for customers after being 
filtered by a rate window of 0.1-yen width. Unfortunately, this raw data is not available to 
us and we cannot use any information regarding the high-frequency raw data. 

Recently, an empirical data analysis by one of the authors of this paper 10| revealed that 
the Weibull probability distribution function, which is defined by 

^m— 1 / 

Pw{t) = m exp I J , (6) 

is a good approximation of the Sony Bank USD/JPY rate in the non-asymptotic regime. We 
should keep in mind that in the asymptotic regime, the FPT distribution obeys the power- 
law — P{t) ~ {a ~ 3.67) — and a crossover takes place at the intermediate time scale 
tc ~ 18,000 [s]. However, the percentage of empirical data for the power-law regime is 
only 0.3%. Similar behaviour was observed for the time-interval distribution between two 
consecutive transactions of Bund futures traded at LIFFE, and the time-interval distribution 
is well-described by the Mittag-LefHer type function [15], namely, the distribution changes 
from the stretched exponential to the power-law for long-time regimes. Therefore, we have 
chosen a Weibull distribution in both asymptotic and non-asymptotic regimes to evaluate the 
average waiting time. The justification for this choice will be discussed later by comparing 
the predicted value with the result from empirical data analysis. 

This Weibull distribution (jH]) is reduced to an exponential distribution for m = 1 and 
a Rayleigh distribution for m = 2. It is easy for us to evaluate the first two moments 
of the distribution Pwii)- We obtained them exphcitly as E(t) = a^/™'r(l + (1/m)) and 
E{t^) = a2/™r(l + (2/m)). Thus, the average waiting time w for the Weibull distribution is 
explicitly given by 

^ = aV™^, (7) 



where T{z) refers to the gamma function, and r{z + 1) = zr{z). In FIG. |3] we plot the 




FIG. 3: Average waiting time w as a function of m for the WeibuU distribution. We set a — 0.5, 1 and 
2. The unit of the vertical axis is second. In the inset, we plot for a — 50.855, which was evaluated from 
empirical data of the Sony Bank USD/JPY rate From the inset, for the parameter m — 0.59, the 



average waiting time w is expected to be w ~ 2534.146 [s] = 42.236 [niin]. This time interval does not differ 
much from the result (w ^ 49 [min]) obtained from our empirical data analysis. 

average waiting time w ([7]) as a function of m for several values of a, and we find that w 
for Sony Bank customers is convenient and reasonable. This is because the empirical data 
analysis reported the value of the parameter as m = 0.59, and the average waiting time w 
for this case is evaluated as w ~ 2534.146 [s] = 42.236 [min] from the inset of this figure. 
From reference information about the scaling parameter a was not available; however, 
we can easily obtain the value as discussed next. 

As mentioned above, several empirical data analyses [loj| revealed that parameter m ~ 
0.59 for the Sony Bank rate. In fact, it is possible for us to estimate a by using the fact that 
the average interval (t) of the rate change is ~ 20 [min ] (TABLE [ID . Next, we obtain the 
following simple relation: 

m = {t). (8) 
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That is {a^/"'/m)T(l/m) = (t). Therefore, 

m{t) 



(9) 



r(-) 

Substituting m ~ 0.59 and (t) ~ 20 x 60 = 1200 [s], we obtain the parameter a for the Sony 
Bank rate as a ~ 50.855. 

It should be noted that the average waiting time is also evaluated by a simple sampling 
from the empirical data of the Sony Bank rates. The first two moments of the first-passage 
time distribution Pwit) — {t) = \imT^^{l/T) Tl=o^ik, (t') = \imT^^{l/T) ELo(^^fc)' 
— are easily calculated from the sampling, and then, the average waiting time is given by 
'it'sampiing = We find that ^sampling — 49 [min] which is not much different from 

that obtained by evaluation {w = 42.236 [min]) by means of the renewal- reward theorem 
([2]). Thus, the renewal- reward theorem introduced here determines the average waiting 
time of the financial system with adequate accuracy, and our assumption of the WeibuU 
distribution for the first-passage time of the Sony Bank rates seems to be reasonable. We 
would particularly like to stress that the 0.3% data from the asymptotic regime does not 
contribute much to the average waiting time. 

A detailed account of this fine agreement between w and ^sampling, more precisely, the 
factors to which we can attribute the small difference between w and ^sampling, will be 



reported in our forthcoming paper 30 1. 



We would now like to discuss the average waiting time for the case in which the first- 
passage process can be regarded as a Poisson process. In the Poisson process, the number 
of events occurring in the interval (0, r] is given by -Pa;(t) = (Ar)'^ e~^'^ /k\ (fc = 0, 1, 2, ■ ■ ■ ). 
For this Poisson process, the time interval t between two arbitrary successive events obeys 
the following exponential distribution: 

P(t) = Xe-^\ (10) 

Then, we obtain the first two moments of the above distribution P{t) as E{t) = 1/A, -E(t^) = 
2/A^. These two moments give the average waiting time for the Poisson process as w = 
E{t^) /2E(t) = 1/A = E(t). Therefore, for the Poisson process, the average waiting time 
w is identical to the average rate interval E{t). This result might be naturally accepted 
because the login-time dependence on the average waiting time is averaged out due to the 
fact that for a Poisson process, each event takes place independently. On the other hand, 



for a non-Poisson process, the login time is essential to determine the average waiting time. 
Thus, the average waiting time becomes E{t^) /2E{t) instead of E{t). 

We should stress that for a given stochastic process whose FPT is a constant T, the 
average waiting time w is easily evaluated hj w = (1/T) sds = T/2. However, the 
result for the above Poisson process implies that for a Poisson process with a mean E{t) 
of the FPT, the average waiting time is not E{t)/2 but E{t). This fact is referred to as 



the inspection paradox 23|] in the field of queueing theory. For a Poisson process, each 
event occurs independently, but the time interval between the events follows an exponential 
distribution. Therefore, there is some opportunity for the customers to log in to their 
computer systems when the exchange rate remains the same for a relatively long time, 
although the short FPT occurs with higher frequency than the long FPT. This fact makes 
w the first moment E{t) of the FPT distribution. 

On the other hand, for the Weibull distribution, the condition under which the average 
waiting time w is the same as the average first-passage time E{t) is given hj w = E{t), i.e. 

2 \ r„ / 1 ^ ^ ^ 



In FIG. m we plot both curves mr(2/m) and {r(l/m)}^, as functions of m. The crossing 
point mx = 1 corresponds to the Weibull parameter for which the average waiting time 
w is identical to the average first-passage time E{t). Obviously, for the parameter regime 
m < nix (the case of Sony Bank), the average waiting time evaluated by the renewal-reward 
theorem becomes longer than the average first-passage time E{t). As mentioned before, 
the Weibull distribution with m = 1 follows the exponential distribution Pw{t) ~ e~*. 
Therefore, as mentioned in the case of the Poisson process, the login-time dependence of the 
customers is averaged out in the calculation of the average waiting time w. As the result, 
w becomes the first moment of the FPT distribution, E{t). However, once the parameter 
of the Weibull distribution has deviated from m = mx = 1, the cumulative distribution is 
no longer an exponential distribution. Thus, the stochastic process is not described by a 
Poisson process and the average waiting time w is not E{t) but w = E{t^) /2E{t) which is 
derived from the renewal-reward theorem. As mentioned above, the empirical data analysis 
revealed that the average waiting time of Sony Bank rates is w sampling = (t^) /2(t) ~ 49 [min]. 
This value is almost twice of the average interval of the rate change data E{t) = (t) ~ 20 
[min] in TABLED This fact is a kind of justification or evidence for the conclusion that the 
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FIG. 4: TOr(2/m) and {r(l/m)}^ as functions of m. The crossing point my, = 1 of the two curves 
corresponds to the WeibuU parameter for which the average waiting time w is identical to the average first- 
passage time E{t). For m < nix (the case of Sony Bank), w > E{t) holds true, whereas for m > rux, 
w < E{t) is satisfied. 

Sony Bank USD/JPY exchange rate obeys non-exponential distribution from the viewpoint 
of the renewal-reward theorem. 

IV. EFFECT OF RATE WINDOW ON FIRST-PASSAGE PROCESSES OF PRICE 
CHANGE 

In the previous sections, we evaluated the average waiting time w for foreign exchange 
rates (price change). In particular, we focused on the average waiting time of the Sony 
Bank rate whose FPT distribution might be represented by a Wcibull distribution. We 
compared the average waiting time w obtained by the renewal-reward theorem with that 
obtained by a simple sampling of the empirical data of the Sony Bank rate and found a 
good agreement between the two. However, as mentioned earlier, the empirical data of Sony 
Bank was obtained by filtering the raw market data using a rate window of 0.1-yen width. 
Unfortunately, the raw data of the real market is not available to us (Sony Bank does not 
keep these records) and the only information we could obtain is the Sony Bank rate itself. 
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Although the main result of the evaluation of the average waiting time is never affected by 
this lack of information, it might be worthwhile to investigate the role of the rate window 
with the assistance of computer simulations. 

For this purpose, we introduce the stochastic process X^g, X^^, ■ ■ ■ , X^^, ■ ■ ■ , where the 
time intervals of two successive time stamps Ati = ki — kQ,At2 = k2 — ki, ■ ■ ■ ,Atn = 
kn — kn-i, ■ ■ ■ obey a WeibuU distribution specified by the two parameters mo and oq. For 
example, for the ordinary Wiener process, the above stochastic process is described by 

Xt+^t = Xt + Yt, Yt=^{0,a^), (12) 
P^{At) = mo ^^^exp fJATl) , (13) 

where Xt denotes the value of price at time t. Obviously, the first-passage time of the 
above stochastic process has a distribution P{t), which is generally different from Pp^(At). 
Nevertheless, here we assume that P{t) also obeys a WeibuU distribution with parameters 
m and a. Next, we investigate the effect of a rate window with a width of e through 
the difference of the parameters, namely, the m-mo and a-ao plots. For simplicity, we set 
a = ao = 1 in our computer simulations. In other words, using computer simulations, we 
clarify the relationship between the input mo and the corresponding output m of the filter 
with a rate window of width e. 



A. The Weibull paper analysis of the FPT distribution 

To determine the output m from the histogram of first-passage time, we use the Weibull 
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28l | for cumulative distribution. Now, we briefly explain the details of the 



paper analysis 
method. 

We first consider the cumulative distribution of the Weibull distributions with a = 1 as 

/oo 
Pwis)ds = exp(-t'"). (14) 

Then, for the histogram of the cumulative distribution obtained by sampling using computer 
simulations, we fit the histogram to the Weibull distribution by using the following Weibull 
paper. 

log{log(i?(t)"^)} = mlogt. (15) 
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Thus, the parameter m of the first-passage time distribution P{t) is obtained as a slope of 
the logt-loglog(i?(t)^^) plot. In the following subsection, we evaluate the parameter m of 
the FPT distribution for the stochastic processes with time intervals of two successive time 
stamps At obeying the Weibull distribution P^(At). Then, we compare mo with m using 
the above Weibull paper analysis (IT^ . 



B. ARCH and GARCH processes as stochastic models for raw market data 

As stochastic model of the raw data of the real market, the Wiener process defined 
by equations (IT^ and (IT^ is one of the simplest candidates. However, numerous studies 
from both empirical and theoretical viewpoints have revealed that the volatility of financial 
data such as the USD/JPY exchange rate is a time-dependent stochastic variable. With 
this fact in mind, we introduce two types of stochastic processes that are characterised by 

25] and GARCH models Q, 27 1 . 



time-dependent volatility — ARCH 



The ARCH(l) model used in this section is described as follows: 

^t+At = ao + alX^, (16) 
Xt+At = Xt+Ar{0,cr^), (17) 



where we assume that the time interval At obeys the Weibull distribution (1131) . From 
the definition in (fT6!) and f|T7|) . we find that the volatility controlling the conditional 
probability density P{Xt) at time t fiuctuates. However, for the limit t —>■ oo, such a 
local time dependence does not prevent the stochastic process from having a well-defined 
asymptotic distribution P{X). In fact, the above ARCH(l) model is characterised by the 
variance cx^ observed over a long time interval t — >■ oo. It is easily shown that cr^ can be 
written in terms of and ai as 

= (18) 

1 — «! 

We choose the parameters ao and ai so as to satisfy a = e = 1. As a possible choice for 
this requirement, we select {ao,ai) = (0.45,0.55). For this parameter choice, the kurtosis 
(Kurt), which is defined by the second and fourth moments of the probability distribution 
function of the stochastic variable X as (X^)/(X^)^, leads to 

6q;? 
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As is well-known, the higher kurtosis values of financial data indicate that the values close 
to the mean and extreme positive and negative outliers appear more frequently than for 
normally distributed variables. In other words, the kurtosis is a measure of the fatness of 
the tails of the distribution. For instance, a normal Gaussian distribution has Kurt = 0, 
whereas distributions with Kurt > are referred to as leptokurtic and have tails fatter than 
a Gaussian distribution. The kurtosis of the ARCH(l) model with {ao,ai) = (0.45,0.55) is 
Kurt = 23. 

Here, we also introduce the GARCH(1,1) model defined as 

c^t+At = ao + aiX^ + pial (20) 
Xt+At = X,+Ar{0,a^), (21) 

where the time interval At is assumed to obey the WeibuU distribution ( fT3ll . The variance 
of the above GARCH(1,1) model observed on a long time interval, t oo, is given by 

i — ai — pi 

To compare the effect of the rate window for the ARCH(l) model with that of the 
GARCH(1,1) model, we choose the parameters ao,«i and (3i so as to satisfy cr = e = 1. 
Among many candidates, we select (ao,ai,/3i) = (0.4,0.3,0.3), which gives the kurtosis 



6aj 

1 - 3al ^2^i/?i - pf 



K^rt = 3 + , = 4.17. (23) 



In FIG. [5l We plot the probability density function (PDF) of the successive increments 
(returns) Y = Xt+At - Xt for the ARCH(l) and the GARCH(1,1) models. The left panel 
is the PDF for the ARCH(l) with (ao, ai) = (0.45, 0.55) and Kurt = 23. The right panel is 
the PDF for the GARCH(1,1) model with (ao,ai, A) = (0.4,0.3,0.3) and Kurt = 4.17. 

For the stochastic processes for modelling the raw real market data, namely, the Winer 
process and the ARCII(l) and GARCII(1,1) models, we determine the parameter m for 
each FPT distribution by means of the WeibuU paper analysis based on f|T5l) . and then 
plot the mo-m relation for each stochastic model. In the left panel of FIG. El we plot 
the iTLQ-m relations for the stochastic processes: the Wiener process and the ARCH(l) and 
GARCH(1,1) models. In this panel, the crossing point between the line m = mo and each 
plot indicates the value of mo for which the distribution P{t) of the first-passage time remains 
the same as that of the time interval At of the raw data P^(At). Therefore, below the line 
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FIG. 5: Probability density function of successive increments (return) Y = Xt+At — Xt- The unit of the 
horizontal axis is yen. The left panel is the PDF for the ARCH(l) model with {aQ,ai) = (0.45,0.55) and 
Kurt = 23. The right panel is the PDF for the GARCH(1,1) model with (ao,ai,/3i) = (0.4,0.3,0.3) and 
Kurt = 4.17. In both panels, the broken line corresponds to the PDF for a normal Gaussian A/'(0, 1). 
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FIG. 6: The mo-m relations for the Wiener process and the ARCH(l) and GARCH(1,1) models (left). 
The right panel shows the mo-dependence of the average waiting time w for each stochastic process. In this 
panel, the unit of the vertical axis is second. Each error-bar is evaluated by 20 independent runs. 



m = tuq, the rate window affects the stochastic process so as to decrease the parameter mo of 
the WeibuU distribution, whereas, above this hne, the parameter mo increases to m (> mo). 
As mentioned earher, empirical data analysis of the Sony Bank rate [ic| suggested that the 
waiting time or the time interval At of the Sony Bank rate obeys the Weibull distribution 
with m = 0.59. Therefore, there is a possibility that the raw data before applying the rate 
window could be modelled by the GARCH(1,1) model. The right panel of FIG. [6] shows the 
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mo-dependence of the average waiting time obtained by the data after the rate window filter 
is apphed. From this panel, we find that the GARCH(1,1) model reproduces the average 
waiting time w for the raw data (WeibuU: exact) more effectively than the other three 
models. 



V. EVALUATION OF EXPECTED REWARDS 

In the previous sections, we considered the average waiting time that a customer has to 
wait until the next update of the Sony Bank rate after she or he logs in to her/his computer 
systems. To evaluate the average waiting time, we used the renewal-reward theorem, which 



is wellknown in the field of queueing theory 23|, |2J]. Besides the average waiting time, we 



need to consider another relevant quantity while investigating the statistical properties of 
the Sony Bank rate from a different viewpoint. For instance, the cumulative return that the 
customers can be expected to obtain during the time interval t, 

N(t) 

R{t) = J2 = - X„, (24) 

n=l 

is one of the candidates for such relevant quantities. In the definition of the cumulative 
return flM|) . N{t) refers to the number of rate changes within time interval t. The return 
Yn is the difference between the rates of two successive time stamps n + 1 and n. Then, the 
long-time average of the cumulative return defined hj R = limt^oo(-R(^)/^) is rewritten as 

R(t) N(t) R(t) 

R = lim = lim ■ (25) 

t^oo t t~.oo t N{t) ^ ' 

Taking into account the following relation : 

lim ^ = E{tY^ (26) 

and the law of large numbers R{t)/N{t) = S{Y), we obtain a long-time average of the 
cumulative return, i.e. the reward rate R as follows: 

R - (27) 

where we define E{t) as the expectation of the first-passage time t and S{Y) as the average 
of the return Y over the probability distribution P{Y). If we set S{Y) = fi and assume that 
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the first-passage time t obeys the Weibull distribution, the above reward rate R is given by 

/im 



R 



(28) 



Obviously, if the distribution of the difference of the rates Y between two arbitrary suc- 
cessive time stamps obeys symmetric distribution around Y = 0, fi becomes zero, and as 
a result, the reward rate R also becomes zero. For performing theoretical evaluations for 
R, namely, to investigate the skewness S'-dependence of the reward rate R, we assume that 
the stochastic variable Y obeys the following skew-normal distribution, wherein the regions 
(—00, — e], [e, 00) are cut off: 



P(F) 



2 

— e 2 

TT 



H{\Y)Q{\Y\-e) 



(29) 



where G(x) is a Heaviside step function and H(x) is defined by H(x) = J^{dt/y/27i) e~*^/^. 
The variable Y cannot take any values lying in the intervals (—00, — e] and [e, 00) because of 
the method of generating of the Sony Bank rate, as we have already explained in the previous 
sections. That is, the Sony Bank rate changes if and only if the difference = — X„| 

becomes larger than the width of the rate window e. In FIG. [TJ we show the skew-normal 
distribution fl29l) for e = — 1 and A = 1. Note that P{Y) becomes a normal distribution for 




FIG. 7: The skew-normal distribution defined by ([29]) with e = 1 and A = — 1. The unit of the horizontal 
axis is yen. 
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the limits A ^ and e — 0. 

For this skew-normal distribution fl29|) . we easily obtain the average = YP{Y) as 



The second and the third moments ^2 = ^^-P(^)) A^s = ^^-P(^) lead to 

/X2 = ^e-#+i7(6), (31) 



^3 = '""""^'^^^'v^ '""""^^^^"'v^Tm^ 




The skewness S* of the distribution P{Y) is written in terms of these moments as follows: 

^ ^ g((y-/i)3) ^ - 3/i/i2 + 2/i=^ 

(j3 CT'^ ' 



where a is the standard deviation a = y?. 

In following, we evaluate the reward rate as a function of the skewness S for the 
parameter values m = 0.59 and a = 50.855. In FIG. [HI we plot them for the cases of e = 0.1 
and 1. From this figure, we find that the reward rate increases dramatically as the skewness 
of the skew-normal distribution (129|) becomes positive. As the skewness S increases, the 
reward rate R saturates at 0.0003 [yen/s] for our model system, wherein the time interval 
of the rate change t obeys the Weibull distribution with (m, a) = (0.59,50.855), and the 
difference of the rate Y follows the skew-normal distribution ( |29l) . In this figure, we also 
find that if we increase the width e of the rate window from e = 0.1 to 1, the reward rate 
R decreases. For companies or internet banks, this kind of information might be useful 
because they can control the reward rate (this takes both positive and negative values) for 
their customers by tuning the width e of the rate window in their computer simulations. 

Moreover, we should note that we can also evaluate the expected reward R = Rw, which 
is the return that the customers can expect to encounter after they log in to their computer 
systems. By combining the result obtained in this section R = 0.0003 [yen/s] and the 
average waiting time w = 42.236 x 60 = 2534.146 [s], we conclude that the expected reward 
should be smaller than R = Rw ~ 0.76 [yen]. This result seems to be important and useful 
for both the customers and the bank, e.g. in setting the transaction cost. Of course, the 
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FIG. 8: The reward rate i? as a function of the skewness S of the skew- normal distribution. The unit of 
the vertical axis is yen per second. We set the parameters m and a for the WeibuU distribution with respect 
to the first-passage time as (m, a) = (0.59, 50.855). We plot it for the cases of e = 0.1 and 1 with A = —1. 

probabilistic model considered here for P{t) or P{Y) is just an example for explaining the 
stochastic process of the real or empirical rate change. Therefore, several modifications are 
needed to conduct a much deeper investigation from theoretical viewpoint. Nevertheless, 
our formulation might be particularly useful in dealing with price changes in real financial 
markets. 

VI. SUMMARY AND DISCUSSION 

In this paper, we introduced the concept of queueing theory to analyse price changes in a 
financial market, for which, we focus on the USD/JPY exchange rate of Sony Bank, which 
is an internet-based bank. Using the renewal-reward theorem and on the assumption that 
the Sony Bank rate is described by a first-passage process whose FPT distribution follows 
a WeibuU distribution, we evaluated the average waiting time that Sony Bank customers 
have to wait until the next rate change after they log in to their computer systems. The 
theoretical prediction and the result from the empirical data analysis are in good agreement 
on the value of the average waiting time. Moreover, our analysis revealed that if we assume 
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that the Sony Bank rate is described by a Poisson arrival process with an exponential FPT 
distribution, the average waiting time predicted by the renewal-reward theorem is half the 
result predicted by the empirical data analysis. This result justifies the non-exponential time 
intervals of the Sony Bank USD/JPY exchange rate. We also evaluated the reward that a 
customer could be expected to encounter by the next price change after they log in to their 
computer systems. We assumed that the return and FPT follow skew-normal and WeibuU 
distributions, respectively, and found that the expected return for Sony Bank customers is 
smaller than ~ 0.76 yen. This kind of information about statistical properties might be 
useful for both the costumers and bank's system engineers. 

As mentioned earlier, in this paper, we applied queueing theoretical analysis to the Sony 
Bank rate, which is generated as a first-passage process with a rate window of width 2e. We 
did not mention the high-frequency raw data underlying behind the Sony Bank rate because 
Sony Bank does not record raw data and the data itself is not available to us. Although our 
results did not suffer from this lack of information, the raw data seems to be attractive and 
interesting as a material for financial analysis. In particular, the effect of the rate window 
on high-frequency raw data should be investigated empirically. However, it is impossible 
for us to conduct such an investigation for the Sony Bank rate. To compensate for this 
lack of information, we performed the GARCH simulations, wherein the duration between 
price changes in the raw data obeys a WeibuU distribution with parameter mo. Next, we 
investigated the effect of the rate window through the first-passage time distribution under 
the assumption that it might follow a Weibull distribution with parameter m (7^ mo). The 
rrio-m plot was thus obtained. Nevertheless, an empirical data analysis to investigate the 
effect might be important. Such a study is beyond the scope of this paper; however, it is 
possible for us to generate the first-passage process from other high-frequency raw data such 
as BTP futures. This analysis is currently under way and the results will be reported in our 



As mentioned in TABLE I, the amount of data per day for the Sony Bank rate is about 
70 points, which is less than the tick-by-tick high-frequency data. This is because the Sony 
Bank rate is generated as a first-passage process of the raw data. This means that the 
number of data is too few to confirm whether time-varying behaviour is actually observed. 
In our investigation of the limited data, we found that the first-passage time distribution 
in a specific time regime (e.g. Monday of each week) obeys a Weibull distribution, but 



forthcoming article 
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the parameter m is slightly different from 0.59. However, this result has not yet confirmed 
because of the low number of data points. Therefore, we used the entire data (about 31, 000 
points from September 2002 to May 2004) to determine the WeibuU distribution. In this 
paper, we focused on the evaluation of the average waiting time and achieved the level of 
accuracy mentioned above. In addition, if Sony Bank records more data in future, we might 
resolve this issue, namely, whether the Sony Bank rate exhibits time-varying behaviour. 
This is an important area for future investigation. 

Although we dealt with the Sony Bank USD/JPY exchange rate in this paper, our ap- 
proach is general and applicable to other stochastic processes in financial markets. We hope 
that it is widely used to evaluate various useful statistics in real markets. 
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